In Search of Microscopic Evidence for Molecular Level Negative Thermal Expansion 

in Fullerenes 
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We report the high-resolution far infrared vibrational properties of Ceo and endohedral Kr@Ceo 
fullerene as a function of temperature. Anomalous softening of the Fi u (l) mode (526 cm -1 ) is 
observed throughout the temperature range of investigation (300 - 10 K) suggesting that the 
fullerene cage may expand at low temperature in these molecular solids. To test this idea, we 
combine these results with temperature and pressure dependent Raman, infrared, and Kr extended 
x-ray absorption fine structure data from the literature to provide a comprehensive view of local 
cage size effects. The results are consistent with a recent molecular dynamics study [Kwon et at, 
Phys. Rev. Lett. 92, 15901 (2004)] that predicts negative thermal expansion in carbon fullerenes. 



PACS numbers: 61.48.+C, 78.30.-j, 65.40.-b, 63.20.-e 
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INTRODUCTION 

The interplay between thermal, mechanical, and vi- 
brational properties that gives rise to negative thermal 
expansion (NTE) in complex materials is an area of sus- 
tained scientific interest. ZrW20s is a prime example, 
where bulk NTE is present over a wide temperature 
range and is attributed to unusual low-frequency dynam- 
ics 0, S 01 • Various layered and molecular solids ex- 
hibit similar behavior although the mechanisms driving 
NTE can be different. For instance, graphite and re- 
lated layered solids display anisotropic negative thermal 
expansion even well above room temperature 0. 0] . In 
graphite, NTE is primarily due to the intra-layer ten- 
sion effect, whereas in molecular solids such as Sm2.72C6o, 
NTE originates from the mixed electronic configuration 
of Sm pj. Based on bonding pattern similarities and 
compressibility arguments, it has been anticipated that 
other nanocarbon materials (including fullerenes, small 
cage compounds, and nanotubes) may also exhibit un- 
usual properties ||. Independent of the many scientific 
discoveries in this area, the technological dream is to ex- 
ploit negative and zero expansion materials in composites 
and molecular electronics devices that have specifically 
tailored thermal expansion coefficients. 

A number of recent molecular dynamics simulations 
have focused on the thermal properties of small cage 
fullerenes and carbon nanotubes d,El>ElE]]- Of partic- 
ular interest is the molecular dynamics investigation by 
Kwon et ai, which predicted NTE in Cgo up to mod- 
erate temperatures due to an unusual combination of 
structural and vibrational entropic effects |8j ■ The calcu- 
lated value of the maximum volumetric thermal expan- 
sion coefficient is-lx!0~ 5 K -1 ||, comparable to that in 



ZrM^Os 0- It is important to note that NTE in Cgo is 
redicted to be a molecular-level rather than bulk effect 
. As a consequence, local or "microscopic" probes are 
required to investigate this prediction. At this time, the- 
oretical results for carbon nanotubes are more controver- 
sial 0, 0, 0, |n| , demonstrating the need for additional 
work in this area. 

At ambient temperature and pressure, Ceo fullerene 
has the shape of a truncated icosahedron (point group I/j) 
with 12 pentagonal and 20 hexagonal faces. It displays 
both long (between a hexagon and a pentagon, 1= 1.45 
A) and short (between hexagons, I— 1.40 A) bonds [l2j |. 
Due to the high point symmetry, isolated Gqq molecules 
display 46 distinct vibrational modes. Of this set, 4Fi u 
are infrared active, and 2A g + 8 H 9 are Raman active. 
These modes can be differentiated by the percentage of 
radial and tangential displacement of molecule in a par- 
ticular vibrational mode [lflj . For example, Fi„(l) and 
Fi u (2) modes are 93.5 and 66.6% radial while Fi u (3) 
and Fi„(4) are 48.7 and 99.5% tangential in character, 
respectively 0, Q. In the solid state, Cqq molecules 
crystallize in face-centered cubic arrangement with space 
group Fm3m {O h 5 ) 15]. The room temperature phase 
is orientationally disordered due to rapid and continu- 
ous rotation of the Cqo icosahedra . Ceo undergoes a 
first-order structural transition to a simple cubic struc- 
ture (Pa3/Th 6 ) at 260 K and to an orientationally or- 
dered phase at around 90 K 0, 0, 0] . Bulk structural 
studies show that lattice parameter trends are standard 
over the full temperature regime (Table [ij jifT ]. 

Incorporation of a noble gas atom such as Kr inside the 
cage leads to a stable endohedral fullerene (Kx@Cp o) with 
an average Kr-C separation of 3.540(3) A jHHg . Com- 
parison of this distance with the sum of the van der Waals 
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TABLE I: 

Structural studies of crystalline Ceo fullerene. 



T(K) 


a (A) 


Ref. 


P (GPa) 


a (A) 


Ref. 


270 


14.1543(2) 


[21] 


0.35(5) 


14.067(2) 


[23] 


170 


14.0708(1) 


f21] 


0.41(5) 


14.031(5) 


f23] 


110 


14.052(5) 


[22] 


0.55(5) 


13.978 


[23] 


11 


14.04(1) 


(19] 


2.4(1) 


13.730 


[23] 


5 


14.0408(1) 


[12 


2.5(1) 


13.633 


[23] 



radii for carbon and krypton (3.72 A) demonstrates that 
the Kr atom is tightly confined within the surrounding 
cage. Recent synchrotron x-ray diffraction work on the 
similar Ar@Cgo indicates that the unit cell parameter is 
slightly larger than that of C« n at 300 K and slightly 
smaller at low temperature j2fj. 

It is clearly of great interest to compare the theoret- 
ical predictions for molecular level NTE in small cage 
fullerenes with experimental observations. In this arti- 
cle, we report a high-resolution far infrared study of Ceo 
and KrQCgo as a function of temperature with particular 
focus on the Fi„(l) and F iu (2) vibrational modes. Un- 
usual softening of the F lu (l) mode is observed with de- 
creasing temperature in both compounds, consistent with 
predictions for an expanded low temperature fullerene 
ball. To further test this idea, we combine our results 
with previous structural and vibrational studies (includ- 
ing variable temperature and pressure Raman, infrared, 
and Kr extended x-ra y a bsorption fine structure data 

m m m m in m in s m m m H3) to evalu- 
ate mode Griineissen parameters and elucidate the ther- 
mal expansion coefficient of C@o- Many of the mode 
Griineissen parameters are negative, and we estimate 
that a V ib is on the order of -10~ 5 K _1 . We also discuss 
the subtle guest-host interaction in Kr@C6o- 
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FIG. 1: (Color online) (a) High resolution far infrared trans- 
mittance spectra of Ceo • Note that the Fi„ (1) mode splits and 
softens at low temperature, whereas the Fi„(2) mode hardens. 
The inset shows a close-up view of Fi„(l) mode softening, (b) 
Close-up view of the low temperature Fi„(l) mode softening 
in Kr@C6o- The inset displays previous Kr-EXAFS results 
(25jl . In all panels, the red (dotted) lines represent 300 K re- 
sults, whereas the blue (solid) lines correspond to 4 K data. 
In the inset, the green (solid) line corresponds to 77 K data. 



METHODS 

Endohcdral fullerenes were synthesized using the high- 
pressure, high-temperature technique described earlier 
|25| . The resulting product contained both Cgo pow- 
der and ~1% Kr@C6o- Multiple separations were carried 
out using a Jasco Gulliver 1500 high-performance liquid 
chromatography system, equipped with recycling, auto- 
matic injection, and a diode array detector (300 - 900 
nm). Using a detailed procedure involving recycling and 
sequential chromatography [25j , about 1 mg of 99% pure 
Kr@C6o was obtained. The pristine Cgo used in this work 
was obtained from Aldrich (99.9% purity). 

Variable temperature far infrared transmittance mea- 
surements were performed using a Bruker 113V Fourier 
transform infrared spectrometer, covering the frequency 
range from 22 - 650 cm -1 with both 0.1 and 0.5 cm -1 



resolution. A helium-cooled bolometer detector was em- 
ployed for added sensitivity. The measurements were 
done using isotropic, pressed pellets (~1% by weight) 
with paraffin as the matrix material. Low temperature 
spectroscopies were carried out with a continuous-flow 
helium cryostat and temperature controller. Standard 
peak fitting techniques were employed, as appropriate. 



RESULTS AND DISCUSSION 

Figure ^shows the far infrared transmittance spectra 
of fullerene C 6 o and Kr@C 60 at 10 and 300 K. For C 60 , 
the radial Fi„(l) and Fi„(2) modes are observed at 526 
and 576 cm -1 , consistent with previous work |38j. The 
most striking result is the temperature dependence of the 
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Fi u (l) mode, which softens by ^0.5 cm" 1 with decreas- 
ing temperature. This effect occurs gradually throughout 
the entire temperature range. Note that the fine struc- 
ture in Fi„(l) is a result of crystal field symmetry break- 
ing below the 260 K transition and was investigated in 
detail by Homes et al. j38[. At low temperature, Fi„ 
splits as Fi„ — > A u + E„ + 3F„, of which only the F u 
modes are infrared active The unusual softening 

of Fi„(l) was not, however, discussed in earlier measure- 
ments. A similar softening effect appears in endohcdral 
Kr@C6o, as shown in Fig. ^b). Taken together, these 
temperature dependent mode softening trends support 
a more relaxed ball at low temperature. In contrast to 
the observed softening in Fi u (l), the other radial mode, 
Fiu (2), exhibits normal behavior, that is, hardening with 
decreasing temperature. The mode sharpens and grows 
in strength as the sample is cooled. The absence of crys- 
tal field- induced fine structure in Fi„(2) has been pre- 
viously attributed to the degree of overlap between 6:6 
bond and the adjacent pentagon 

Figure displays the peak position of the Fi„(l) mode 
in pristine Ceo and Kr@Ceo as a function of temperature. 
Two major trends are observed. First, the Fi u (l) mode 
softens throughout the measured temperature range for 
both Ceo and Kr@C6o- Second, the Fi u (l) mode hard- 
ens by ^2 cm" 1 upon endohedral incorporation 0]. A 
hardening trend was also observed for the Fi„(l) mode 
in a previous study [i^. The Fi„(2) mode frequency in 
Kr@C6o (577 cm -1 ) is also ~1 cm" 1 higher than that 
in the pristine material. Noble gas-containing endohc- 
drals are traditionally considered to be van der Waals 
complexes due to the inert nature of noble gas guest 
atom. In such a picture, only modest spectral changes 
(such as the appearance of a rattling mode) are antici- 
pated upon endohedral incorporation. The observed ra- 
dial mode hardening (Fig. |2J) and tangential mode soften- 
ing 0] in Kr@C6o compared with pristine Cgo suggests 
a modest guest-host interaction and a slight modification 
of the encapsulant/cage interaction potential. 

Returning our attention to the temperature dependent 
mode softening problem, a literature survey shows that a 
number of Raman active modes also display unexpected 
low temperature softening [23, US E p .l 3 (M . Just to cite a 
few examples, van Loosdrecht et al. |2JJ observed soften- 
ing of several different modes (H g (3), H ff (4), A g (2) and 
H ff (8)) both above and below the 260 K transition tem- 
perature. (All modes harden at the transition temper- 
ature.) Hamanaka et al. showed that the H s (7) mode 
softens below 260 K, whereas the A g (2) mode softens be- 
low 160 K. The H s (3) and H g (5) modes remain almost 
unchanged with decreasing temperature [2^]. In contrast, 
Horoyski et al. reported almost no change in A g (2) and 
H g (4) in the lower temperature range 01 ■ Dong et 
al. observe higher order overtone and combination modes 
in Cqo, although their exact temperature dependence is 
not clear [4l|. Taken together, we find that several (but 
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FIG. 2: (Color online) Temperature dependence of the center 
peak position of the Fi„(l) mode in Ceo and Kr@Cso. The 
blue-shift of this mode in the endohedral material suggests 
modest guest-host interactions. Two independent runs are 
shown for Ceo- Error bars are indicated. 



not all) Raman investigations support unusual low tem- 
perature vibrational mode softening. 

It is of interest to understand the microscopic origin of 
the unusual mode softening and its implications for the 
bulk properties. In general, mode softening can be cor- 
related to two major effects: volume and charge. The 
charge effect is well-known in organic charge transfer 
salts ^3 EH ■ Typically, volume contributions dominate 
in high pressure experiments. High-pressure vibrational 
studieSjdetailed bel ow also indicate extensive mode soft- 
ening [Hjm000IH,li3- Hence, we attribute both 
pressure and temperature dependent mode softening to 
cage volume effects . 

Volume change can be quantified by the thermal ex- 
pansion coefficient as 



lavCy 

V m B 



XT \ - 

v m 



(1) 



Here, B (= 1/xt) is the bulk modulus, ^ av (= 
Si I Tlii °i) i s the average of mode Griineisen parame- 
ters 7i weighted by the mode specific heat a , C'v is total 
molar specific heat at constant volume, and V m is the 
molar volume |5Cj . The mode Griineisen parameters are 
dimensionless quantities that relate the fractional change 
in volume to the fractional change in frequency of a given 
mode. They are defined as 



7» = 



/ 6 In Ui \ 




' '8u>i\ 


\5lnV J 







(2) 



where u>i is the frequency of the i mode. At high 
temperatures, the j av can be taken as the arithmetic 
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mean of mode Griineisen parameters 7*. Note that ther- 
mal expansion will be positive or negative depending 
upon whether positive or negative 7^ predominate in the 
weighted average [jl lil loOf . 

To evaluate these quantities, we examined the pres- 
sure dependent vibrational properties data that are avail- 
able in the literature. Meletov et al. measured the Ra- 
man spectrum of Ceo between 200 - 800 and 1350 - 1700 
cm -1 [3lLl51| . The results clearly show three well-defined 
phases: below 0.4 GPa, between 0.4 and 2.4 GPa, and 
above 2.4 GPa. These pressure-induced phase transi- 
tions are attributed to orientational ordering from face- 
centered cubic to a simple cubic structure and finally to 
a rotation-free orientationally-ordered simple cubic phase 
in the high pressure regime pilj . Comparison of lattice 
parameter trends with temperature and pressure (Ta- 
ble P) demonstrates the rationale for concentrating our 
analysis on the pressure range below 0.4 GPa, although 
to obtain the fullest possible mode contributions, we also 
employ higher pressure results in our secondary analysis, 
below. Returning to the data of Meletov et al. , al- 
most all the modes soften up to 0.4 GPa, whereas H ff (3), 
H s (4) soften through the higher pressure range as well. 
Similar softening has been observed for the A s (2) mode 
at 0.35 GPa |22j. Further, Snoke et al. showed that 
the Fi u (l), H ff (3), and H 9 (4) modes soften up to 14 GPa 
[3^ |. Unusual softening of Fi„(l) and several Raman and 
infrared inactive modes are also observed in other inde- 
pendent variable pressure measurements 000113 
Martin et al. highlight pressure-induced softening of 
many inactive modes up to 2.5 GPa as well 0]. Despite 
these experimental findings, no comprehensive explana- 
tion for the unusual mode softening of Fi u (l) and the 
aforementioned Raman active modes has been advanced 

m 

Table[n]summarizes our extended analysis of the afore- 
mentioned pressure dependent vibrational studies of Cgo, 
undertaken as part of our search for microscopic evidence 
of molecular level NTE in fuller ones. As discussed pre- 
viously, many modes soften with increasing pressure (es- 
pecially between - 0.4 GPa), giving negative Griineisen 
parameters for the majority of modes. The size and sign 
of the Griineisen parameters are compatible with small 
but significant expansion in intramolecular bond len gth s 
as pressure is increased or temperature is decreased |49J . 
Also, note that both radial and tangential modes exhibit 
pressure- induced softening [55| . The average value of the 
mode Griineisen parameters including most of optically 
active modes (upper portion, Table Hl)l is negative in this 
pressure range (j a v ~ -0.19). In order to further pin 
point the total average mode Griineisen parameter for 
Ceo, we extended this analysis to include several optically 
inactive and two Raman active modes over a somewhat 
larger pressure range 0.5 - 2.5 GPa 35|. The selected 
modes are shown in the lower portion of Table ITT1 Thus, 
combining the data for these modes in the higher pres- 



TABLE II: 

Mode Griineisen parameters of the Raman, infrared active, 
and inactive modes of C60 in two different pressure ranges: 
- 0.4 and and augmented to include selected data from 0.5 - 

2.5 GPa. The reported value of xt = 6.9 x 10" 2 GPa -1 is 
used in the calculation |H3| . 

Mode uj doj/dP Ref. ji 
(cm- 1 ) (cm" 1 /GPa) 



H 9 (l) 


272 


-15.00 


f31] 


-0.80 


H 9 (2) 


435 


-12.00 


f31] 


-0.40 


A 9 (l) 


496 


-12.00 


[31] 


-0.35 


Fm(l) 


526 


-1.60 


f37] 


-0.04 


Fi„(2) 


576 


2.40 


[37] 


0.06 


H 9 (3) 


710 


-12.00 


[31J 


-0.24 


H 9 (5) 


1099 




[31]" 




Fm(3) 


1182 


3.90 


[37] 


0.05 


H 9 (6) 


1248 




[31]° 




Fi« (4) 


1429 


6.10 


[37] 


0.06 


A 9 (2) 


1467 


-18.00 


[31] 


-0.18 


H 9 (8) 


1570 


-7.00 


[31] 


-0.06 



7a „ ~ -0.19 6 



H 9 (4) 


772 


-2.70 


[36J C 


-0.05 


H 9 (7) 


1422 


9.80 


[36] r 


0.10 


H„(3) 


668 


0.62 


[35] c 


0.01 


F 2u (2) 


712 


-1.70 


[35] c 


-0.04 


G„(2) 


741 


-1.90 


[35] c 


-0.04 


F 2u (3) 


796 


0.73 


[35] c 


0.01 


H„(4) 


824 


3.00 


[35] c 


0.05 


G 9 (2) 


961 


1.80 


[35] c 


0.03 



Jav —0.W d 



"Ref. 121 

'Low pressure regime: - 0.4 GPa. 
c Slope between 0.4 - 2.4 GPa 

^Inclusion of Raman and infrared active modes from the low 
pressure regime as well as selected inactive modes with pressure- 
dependent frequency shifts obtained from the high pressure regime, 
as indicated in the Table. 



sure range (0.5 - 2.5 GPa) with the aforementioned mode 
Griineisen parameter of the optically active modes in the 
lower pressure range yields an averaged value of the mode 
Griineisen parameter that is also negative ( , y av ~ -0.10). 

Can we use this information to calculate the thermal 
expansion coefficient, a? Certainly, the negative value of 
the averaged Griineisen parameter is directly related to 
a negative thermal expansion coefficient. With precise 
knowledge of mode Griineisen parameters and phonon 
densities of states for each mode, the calculation of ther- 
mal expansion coefficient is straightforward (Eqn. (1)). 
It is, however, important to point out that isolated Cgo 
has 46 vibrational modes and that we are able to evalu- 
ate Griineisen parameters for only a small subset based 
upon currently available data, despite the two different 
trials presented in Table [H] To finesse this problem, we 
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designate a V ib as the vibrational part of the thermal ex- 
pansion coefficient for the optically-active modes in the 
low pressure range (0 - 0.4 GPa) shown in Table ITU Us- 
ing the reported density of states HfHSlJj each mode was 
weighted to calculate the total specific heat from mode 
specific heat and finally our pseudo thermal expansion 
coefficient, a V ib 0]. Based upon this data, we find a v u, 
to be on the order of -1CP 5 K _1 . Inclusion of the inactive 
modes over the larger pressure range (bottom portion of 
Table ITTjl does not modify this general picture. Note that 
this value quantifies molecular rather than conventional 
lattice effects. 



Independent experimental support for this vibrational 
analysis and the predictions of Kwon et al. comes 
from our reanalysis of the extended x-ray absorption fine 
structure study of Kr@C6o, performed at the Kr-edge of 
the encapsulated Kr atom at 77 and 300 K j^. X-ray 
absorption spectroscopy is a sensitive and direct probe 
of molecular size and other local effects. The inset of 
Fig. njb) shows the Kaiser-Bessel Fourier transformed 
data. A small shift in the peak position (to higher values 
of R) is clearly evident at low temperature and is con- 
sistent with the mode softening of various infrared and 
Raman active modes. The Fourier analysis and fitting 
is tricky, but the authors extract a Kr-C path distance 
of 3.540(7) A at low temperature and 3.537(10) A at 
300 K, indicating the possibility of a slightly relaxed low 
temperature cage Translating this observation into 
a thermal expansion coefficient, we obtain on the order of 
-10~ 5 K" 1 , similar to that extracted from our vibrational 
analysis. 



Combined, these experiments provide broad support 
for the proposed molecular-level NTE in small carbon 
cage compounds. In Cgo and Kr@C6o, selected modes 
soften as temperature is lowered and/or pressure is in- 
creased. The mode softening correlates directly with 
the size and geometry of the fullerene cage and NTE 
at molecular level offers a universal explanation of both 
temperature and pressure effects 0, |55|. The trend 
holds even in the presence of subtle guest-host inter- 
actions in Kr@C6o- As discussed by Kwon et al. @, 
low-temperature cage expansion may originate from an 
unusual combination of structural and vibrational en- 
tropic effects. Hence, NTE may be an intrinsic prop- 
erty of fullerene and other small-cage molecules, although 
in bulk measurements, the molecular-level effect may be 
dominated by bulk lattice expansion. On the other hand, 
molecular electronics applications exploit single molecu- 
lar properties. Understanding the thermal properties of 
functional, assembled molecular-based materials presents 
many future opportunities and challenges. 



CONCLUSION 

We report the the high-resolution far infrared vibra- 
tional properties of Ceo and KrQCgo- We observe that 
the Fi M (l) mode softens with decreasing temperature 
in both materials (300 - 10 K) and propose that the 
softening is correlated to low temperature expansion of 
fullerene cage. This finding is consistent with previ- 
ous temperature/pressure dependent vibrational spec- 
troscopy and is also supported by our reanalysis of tem- 
perature dependent Kr extended x-ray absorption fine 
structure data jH E3, El S El El iHl HI iHl El IS 
EJ. To further test this hypothesis, high pressure spec- 
tral results 0, El EJ were employed to calculate the 
mode Gruneisen parameters. Many of these Griineisen 
parameters are negative. These results are consistent 
with the recent prediction for molecular-level NTE by 
Kwon et al. Further work is needed to extend these 
ideas to doped fullerenes and other curved molecular sys- 
tems. 
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